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The anisotropic Kondo necklace model in 2D and 3D is treated as a genuine model for magnetic 
to Kondo singlet quantum phase transitions in the heavy fermion (HF) compounds. The variation 
of the quantum criticalpoint (QCP) with anisotropy parameters has been investigated previously 
in the zero field case [l|. Here we extend the treatment to finite fields using a generalised bond 
operator representation including all triplet states. The variation of critical t c with magnetic field 
and the associated phase diagram is derived. The influence of anisotropies and the different g- 
factors for localised and itinerant spins on the field dependence of t c is also investigated. It is 
found that three different types of behaviour may appear: (i) Destruction of antiferromangetism 
and appearance of a singlet state above a critical field, (ii) The inverse behaviour, namely field 
induced antiferromagnetism out of the Kondo singlet phase, (iii) Reentrance behaviour of the 
Kondo singlet phase as function of field strength. 

PACS numbers: 75.10 Jm, 75.30 Mb 

I. INTRODUCTION 

The Kondo necklace (KN) type models are useful to discuss the quantum phase transitions between Kondo singlet 
and antiferromagnetically ordered states such as found in heavy fermion compounds 0, H, EL [H, [ll- They have 
been originally proposed by Doniach Q for the one-dimensional case as a simplified version of the itinerant Kondo 
lattice (KL) models Q. Thereby the kinetic energy of conduction electrons is replaced by an intersite exchange 
term. For a pure xy-type intersite exchange this may be obtained by a Jordan- Wigner transformation. However in 
higher D the replacement cannot be justified easily. The intuitive argument is that at low temperatures the charge 
fluctuations in the Kondo lattice model are frozen out and the remaining spin fluctuation spectrum can be simulated 
by an antiferromagnetic inter-site interaction term of immobile r spins coupled by a Kondo interaction to the local 
noninteracting spins S. Recent exact diagonalisation studies on finite clusters for both Kondo lattice and Kondo 
necklace models have indeed found that the competition between on-site Kondo singlet formation and AF inter-site 
correlations are very similar in both models [3] • A more formal way to get rid of charge fluctuations in the KL 
model is an inclusion of a Hubbard term for conduction electrons (KLU model) which leads to the isotropic KN model 
in the large U limit at half filling. 

Nevertheless one should consider the Kondo necklace model for D^2 as a model in its own right which is suitable 
for studying quantum phase transitions between a Kondo singlet (KS) and an antiferromagnetic (AF) phase. In its 
original form the local Kondo exchange is isotropic while the intersite exchange is of xy-type. This model has U(l) 
symmetry. Later more general models with arbitrary anisotropy of Kondo as well as intersite exchange terms have 
been considered [ll[ . Indeed compounds which exhibit the KS to AF quantum phase transition have mostly uniaxial 
symmetries. A full account of the influence of uniaxial anisotropies of both terms in the Kondo necklace Hamiltonian 
on the quantum critical point has been given in Ref. [l|. The D=2 KN model without any anisotropy may also be 
understood as a special case of a bilayer Heisenberg model [T3, EH where the intersite bonds are cut in one of the 
layers. A reintroduction of holes in this case leads to the KNtJ model which is related to the KLU model away from 
half filling [HEl- 

In the general anisotropic KN model the quantum phase transition is achieved by varying the ratio of 'hopping' 
t, i.e., the intersite interactions of r spins to the on-site Kondo interaction J. In practice this is achieved by varying 
pressure (hydrostatic or chemical). An alternative way to arrive at the QCP is to apply an external magnetic field 
which breaks the local Kondo singlets and leads to a field dependence of the critical t c . Starting from a noncritical or 
above critical t at zero field the system may then be tuned to the QCP by varying the field strength. This is indeed a 
practical method frequently applied [To| . To investigate field-induced quantum phase transitions from Kondo singlet 
to AF ordered state or vice versa we have extended our previous work to include the effect of the magnetic field 
within the Kondo necklace model. However, this introduces an additional parameter, namely the ratio of effective 
g- factors for the local Kondo spins and the interacting spins. They can be different due to the different strength of 
spin-orbit coupling and crystalline electric field effects involved in the formation of the r and S (pseudo-) spins. A 
mainly numerical study of the itinerant isotropic KL model in a magnetic field has been given previously [16| . 

In Sec. [TTJ we define the anisotropic KN model in a magnetic field and in Sec. IIIII briefly discuss the local state 
space of the Kondo term in an external field. In Sec. II VI we perform the transformation from spin to bosonic variables 
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and in Sec. [V] derive the selfconsistent equations for the mean field boson condensate amplitudes in the singlet and 
antifcrromagnetic phases.. The influence of higher order terms in the Hamiltonian is discussed in Sec. IVII In Sec. IVIII 
we investigate the numerical solutions and discuss the resulting quantum critical t-h phase diagram. A discussion is 
provided in Sec. I Villi and Sec. IIXI finally gives the conclusions. 



II. ANISOTROPIC KONDO NECKLACE MODEL IN AN EXTERNAL FIELD 



To investigate field induced quantum critical behaviour we start from the anisotropic Kondo necklace (KN) model 
with U(l) symmetry [l[ where the field is applied along the anisotropy (z) direction: 

H = Ht + Hj + Hz 

= t J2 + + + ^E« S « + + A<^) + 7 ^ 5^ + 7 ' ^ (1) 

<n,m> n i i 

Here, the summation over nearest neighbor (n.n.) is indicated by brackets and r" is the a-component (a = x,y,z) 
of the 'itinerant' electron spin at site n whereas S"" is the a-component of localized spins at position n. For the 
exchange coupling between the itinerant and localized spins we generally use J x = J as reference energy scale in all 
figures except when stated otherwise. The local anisotropy parameter A is defined by the relation J z =/S.i x between 
the z-axis and in-plane (xy) local exchange. The hopping parameter of the itinerant electrons is proportional to t with 
the anisotropy in the z-direction given by 5. The present model has three control parameters: t/J x and the anisotropy 
parameters (S,A). In the Zeeman term we defined 7 = —g s h and 7' = —g T h with h = /x^H where H is the strength 
of the applied field. Furthermore g s and g T are the gyromagnetic ratios for localised (S°) and itinerant (r") spins 
respectively. They are determined by the combined effect of spin-orbit coupling and crystalline electric fields which 
depends on the degree of localisation or itineracy of electrons. Therefore g s and g r in general need not be equal. We 
shall consider two extreme cases, namely (g s , g T )=(2fl) and (g s , g T )=(2,2). The former is more realistic, since in real 
heavy fermion compounds most of the magnetic response is due to the localised electrons with pseudo-spin S. As in 
previous work p], [T3| the present study of field induced quantum phase transitions in the anisotropic KN model is 
based on the bond operator formulation. Its Hilbert space is spanned by local singlet-triplet states of (S n ,r n ) spin 
dimers represented by bosonic degrees of freedom. In applying this technique to the finite field case we will largely 
use the same or similar notations as in the previous zero field case [1] for consistency. 



III. LOCAL LEVELS AND ZEEMAN SPLITTING 



Before we perform the transformation to bosonic variables it is useful to have a clear understanding of the local 
singlet-'triplet' level scheme as function of anisotropy and magnetic field because certain tendencies of the quantum 
critical behaviour are correlated with the splitting of the ground state and the first excited state of a local bond. 
Therefore we first diagonalise the local Hamiltonian Hl — T~tj + 'Hz- This leads to eigenvalues e i =F li /( J x /4) (i=l-4) 
given by 

ei, 2 = A±7 + , 

e 3 .4 = ±(4 + 7^)5 - A, 

% = ~jm 4^^= GOT' (2) 

with 7± = ^(7 ± 7') = —g±h and g± = ^(g s ±g T ). Defining S* —t z +S z then ei.2 are the energies of the triplet states 
I TT))| 14) with S* = ± 1 respectively which exhibit the linear Zeeman effect. Furthermore 63^4 correspond to the 
triplet Si = state ^(| fl) + I -IT)) an( l singlet state -^(\ |4) — | J,|)) respectively. They show level repulsion in a 

magnetic field, except for g s —g T when the total S* -component is conserved and commutes with TLl- It is obvious that 
the spectrum of eigenstates in Eq.([2]) is invariant under the transformation (g s ,g T ) — > (g r ,g s ). The variation of zero- 
field energy levels with A is shown in Fig. Q] (top) from the xy-U(l) limit (A=0) via the Heisenberg SU(2) point(A=l) 
to the Ising Z2 limit (A=oo). The field dependence for a few selected A-values and g-factors ( g s ,g r ) = (2,0) are shown 
in Fig. [1] (bottom). For ( g s ,g T ) = (2,2) (not shown in Fig. [T]) the ground state and first excited state levels cross at 
h C r/Ja:=0.25, 0.50 for (5, A)=(0,0) and (1,1) respectively. In such a case one may expect that the quantum critical 
t c (h cr .) for the transition from Kondo singlet to AF phase vanishes. 
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FIG. 1: Top: Dependence of local energy levels after Eq. |(2]) on the local anisotropy A (xy-side) or A'=l-1/A (Ising side) at 
zero field. Numbers in parentheses denote the degeneracy of each level. Bottom: Field dependence (h=/i_eH) of local energy 
level for three extreme cases. Note that in the xy-case the excitation energy from ground state to first excited state vanishes 
asymptotically (h^Jz) whereas in the other cases it reaches a constant (HAF) or is equal to a constant (Ising). Here we used 

(g.,gr) = (2,0). 



IV. DERIVATION OF EFFECTIVE BOSONIC MODEL 

We apply the transformation from spin variables (S a ,T a ) to bond variables (s,t a ) (a =x,y,z) to the Hamiltonian in 
Eq. Q. It is given by [Hj]: 

where e a f3j is the fully antisymmetric tensor. By construction the singlet (s) and triplet (t a ) operators generate the 
local eigenstates of Hl for zero field and no anisotropy : 

s f |0) = ^=(| U) - 1 IT)); 410) = -^(1 TT> - I U», 

410) = ^=(1 II) + 1 IT)); 41°) = ^(1 tt) + 1 ID). (4) 
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Here the first and second arrows indicate the z-component of t and S spins respectively. The singlet and triplet 
operators satisfy the usual bosonic commutation relations according to 

[s n , 4] = 1; [W, tl >/3 ] = 6 a ,p; [s n , tl i0 ] = 0. (5) 

All other commutators vanish. The physical states have to satisfy the local constraint s^s n + Y^a Ai a^n,a — 1- This 
transformation leads to an effective Hamiltonian in terms of singlet (s) and triplet (t a ) bosons. It has been argued in 
[H Il7j that one may restrict to the terms which are bilinear in the triplet bosons and we will later in Sect . IVT1 discuss 
to which extent this is justified. For the moment we restrict to the bilinear contribution. It can be written as the sum 
of a local term TLl and an inter-site interaction term TL\. Here 

H L =Hj + H z = y ( - ( 2 + A ) s » s » + ( 2 - A ) f « ,**».« + A «xW + tl iy t n ,„) 

n 

+7— ^ ^ ( S n^ nz tnz s n) ~ *7+ ^('L'«!/ — ^liy^nx) i (6) 
n n 

contains the on-site Kondo interaction {TLj) and the Zccman term (Tiz)- The bilinear interaction term is given by 

<n,m> a—x,y 

+ ^ ] ^ s n^n,z {Smtm.z + ^m.z s m) + ft-C.J . (7) 

<n,m> 

The physical constraint on the local Hilbert space is enforced by adding a Lagrange term at each site with an associated 
chemical potential /j, n , leading to 

H = H L + J2»n(sis n + J2 - !) + M 2) =n Q +nf\ (8) 

n a—x,y^z 

This Hamiltonian is diagonalised within a mean field approximation for the bond-operator singlet and triplet bosons. 
We assume that in general there are three independent bosonic amplitudes which characterise the phases in the (t,h)- 
plane: Singlet s — (s n ) denoting the strength of local singlet formation, staggered triplet t = ±(t n>x ) which determines 
the AF order parameter (polarised along x) and homogeneous triplet to = (t n>z ) which determines the magnetisation 
caused by the external field along z-direction. For technical reasons it is of advantage to transform to circular polarised 
transverse triplet coordinates (u n ,d„) with respect to the field direction (z axis). The transformation and its inverse 
are given by 



1 1 



1 i 

In circular triplet coordinates the Hamiltonian may then be written as 

n ° = E{[-t( 2 + A )+^] s »s„ + [ : |(2- A) + f i n ]ti z t nz + 

n 

[^ A + "f++ VnW n U n + [-?-A - 7+ + Vn]d) n d n + 7-(4*rw + *L S «) ~ Mn } , (10) 



and 



= -- [(s^s^Undm + d n u m ) + h.c.) - (s] l s m (u„u, t rl + d n d] n ) + /i.e.)] 
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(n,m) 

t 



+ T [(sUlntnztmz + h.C.) + (s+SmWL + h.C.)]. (11) 



(n,m) 
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Separating the mean values and the corresponding fluctuations the Fourier components of singlet and triplet operators 
are then given by 



s k = VNs, 

u k = VNuS kiQ + ii k , 

d k = ^d5 k , Q +d k , (12) 
tk.z = vNioSkfi +h,z- 

Here we assumed that if AF order appears it will be of the in-plane type (_L z,H). Therefore the case of AF order 
will only be considered for <5 < 1. Instead of to it will later be convenient to use fo which is defined by the relation 
to = s?o- In the paramagnetic phase at zero field only s will be different from zero with a value close to one 0. As 
the field increases s will decrease and simultaneously the triplet amplitude to associated with the uniform induced 
moment increases. If an AF transition occurs, s will further decrease and u,d will increase accordingly with the 
field. The constraint s 2 +t^ + (\u\ 2 + \d\ 2 ) ~ 1 is always respected for any field strength. In the limit h/J x ^S>1 the 
amplitudes s and to become asymptotically equal. On the mean field level this signifies that the ground state of the 
local level scheme becomes an equal amplitude mixture of the two states with S*=0 which show the level repulsion 
in Fig. [T] Therefore it is possible to use the zero-field singlet-triplet bosons in Eq. ^ as a basis even in the present 
finite field problem. On the mean field level, the physical constraint automatically takes care of the change in the 
local ground state wave function. Inserting the above expressions into the Hamiltonian of Eqs. (| 1011 1 [) one obtains a 
bilinear form in the triplet fluctuation operators (u k , d k , t kz ) which may be diagonalised with two separate Bogoliubov 
transformations given for the z-polarisation by 

a k = cosh(6 ktZ )i kyZ + sinh(0 fciZ )it fc z , 



a 



smh(6 kiZ )t k ,z + cosh(#fc iZ )f_j, z , (13) 



and for the two circular polarised triplets by 



A k = cosh(6 k , ± )u k + sinh(6> fei _L)rfi fe , 
B\ = sinh(9 k ,±)uk + cosh(6> fci _i_)cf 



k ' 



B k = cosh(0 fe ,_i_)4 +sinh(0 fe ,_ L )?il i ., (14) 
A^_ k = sinh(0 fc .jj(4 + cosh(0 fej _ijw^ fc . 



The transformation angles k ^ z ,8 kj ± are obtained from the diagonalisation conditions 



tanh20 fc ^ = ; tanh20 fe a = — . (15) 



d z {k) ' d±(k) 
where the longitudinal auxiliaryy functions i z and d z arc defined by 



t s 2 2 7 — 7 7s 2 

f z (k) = —S 7 (k); d z {k) = y.+ X A z +— 6 7 (k). (16) 



and for the in-plane case the auxiliary functions fj_ and dj_ are given by 



.,2 



f±(k) = \l{k)-, d Utd (k) = /i + ^ + \l(k) ± 7 +, 

d±(k) = ^+rfd)=M+j + y7W. (17) 

Here j(k) = X^iLi cos(fcj) with 7(0) = | denotes the n.n. structure factor in dimension D = | (z = coordination 
number of the simple cubic lattice). It should not be confused with the Zccman energies r y± defined in Sect. IIIIl The 
Bogoliubov transformations in Eq. (|13p and Eq. ([HI) yield the diagonalised bilinear Hamiltonian 

n mf = Eo + 5}n A (fc)At(*)A(A) + n B {k)B\k)B(k) + n z {k)a\ ak \. (is) 
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where the triplet mode frequencies f2 Q (fc) (a=A,B,z) in the mean field Hamiltonian are given by 

Qa (k) = uJ A (k) +7+, 
il B (k) = uJ B (k) - 7+, 

Q z (k) = LJ z {k), (19) 

UJ A {k)=LO B {k) = [d ± (kf -4.f ± (fc) 2 ]5 =LJ ± (k), 

io z {k) = [d z {kf -4f z (k) 2 }^. 

The excitation energies f2 Q (fc) depend on the field both explicitly through 7+ (for a=A,B) and implicitly via ui a (k) 
which is determined by the field dependent singlet and triplet amplitudes. In the nonmagnetic phase (u, d, io=0) the 
singlet-triplet excitation gap is given by the minimum excitation energy at the incipient ordering wave vector. For 
D=2 this is at Q=(7r, tt), explicitly E g — min{tt a (Q), a — A,B,z}. In the approach from the nonmagnetic side the 
quantum critical line t c (h) is then defined by the vanishing of E 5 . 

The ground state energy Eo is a function of three control parameters (t/J x ,A,d), four singlet-triplet expectation 
values (s, to } u, d) and the chemical potential [i. Writing the transverse mean values explicitly in terms of amplitudes 
and phases according to 

U = uexp(i$ u ), d — dexp(i^d), (20) 
the ground state energy can be written as 

E (^-,A,5;s,u,d,i Q ) = N[-\{2J X + J z )si 2 + [is 2 - n + 

J 1 Z 

(-j + A - -^zts 2 )(u 2 + d 2 ) + 7+(u 2 - d 2 ) + -ts 2 ud cos($ M + $ d ) (21) 
+ (i(2J, - J z ) + fj, + ^zSts 2 )P + 2 7 _s< ] 



2 

ka 

As a first step we determine the triplet condensate phases (§ M ,®d) by minimization. Since only one term depends 
on the phases we obtain the extremal condition sin(<I> ll + <I>d)=0. The minimum must also satisfy cos($ u + &d) < 0. 
This is achieved for <j> u + = mr with n an odd integer, i.e. the sum is only determined modulo 2tt. Since there 
is no condition for the difference of phases one of them is arbitrary. We then choose ($ u , <J>d)=(0,7r). The remaining 
continuous degeneracy with respect to the phase difference is a signature of the Goldstone mode which is present 
throughout the AF phase. 



SELFCONSISTENT EQUATIONS FOR SINGLET AND TRIPLET AMPLITUDES AND 

MAGNETISATION 



The minimisation of the ground state energy Eo in Eq. (|21[) leads to selfconsistent coupled equations for the 
condensate amplitudes s, u, d, to and the chemical potential \x. Their structure is slightly different in the nonmagnetic 
(t=0) and magnetic (i > 0) case, therefore we write them explicitly for both. For a convenient expression of the 
selfconsistency equations we define the quantity fo by the relation to = fos. Furthermore we introduce the Brillouin 
zone integrals F a and G Q (a=A,B,z) given by 

F = — da 

A ~ J U>ka ' 

K 

t X - dq(fc)-2/ a (fc) _ -a 

N ^ LOka 

K 

which appear in the extremal conditions when differentiating the last term in E with respect to /1 or s respectively. 
Here we have defined 7^ = 7^ S = 7fc and 7| = 6-fk- Note that even in the case of finite field (7± 7^ 0) it is the 
frequency LOka and not the mode frequency Qka which appears in the expressions for F Q and G Q . 
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FIG. 2: Singlet (s) and triplet (io,t) amplitudes and their associated uniform (mo = sio) and staggered (m s — st) moments 
as function of field strength normalised to the critical field in the KS and AF region. Here h c (S, A)/ J x ~0-55. At the critical 
field (io/s) = 0.25. This case with t/t c (<5, A)=0.6 (t/J x =0.517) corresponds to the the upper curve of Fig. [4] (top). Here 
( 5s)5t )=(2-0). 



Nonmagnetic Kondo spin singlet phase in external field 

In this case the transverse triplet amplitudes vanish, i.e., u=d=0 and the minimization of Eo with respect to 
(s, to = sfo,/x) leads to the selfconsistent set of equations 



.s 2 



1(5 -E^a) 



f = ^ / 23 ^ 

2ii+\{2J x - J z )+z5ts 2: 1 ' 



Here s, /i are found by iteration which determines tq and hence to completely via the second equation above. For zero 
field, i.e. 7_ = 0, the induced fo vanishes and then Eq. reduces to the paramagnetic case of Ref.[l|. 

Antiferromagnetic phase in external field 

The additional minimization with respect to u and d amplitudes leads to a direct relation for the chemical potential 
fi in terms of the singlet amplitude s which is a generalised version of the zero-field AF case [l| : 

tih) = \ztl?-± + [( 1 - z tsr+^- (24) 

At this stage it is convenient to introduce an auxiliary parameter k which controls the effect of the external field on 
the boson amplitudes; it is defined by 



-WHl + bg^ + Tg^ (25) 

Obviously k(0)=1 and n(h) < 1 since 7+ < 0. Using the expression for the zero-field chemical potential /Lto=/i(0) we 
may also write 

z 9 (k — l) 2 

H{h) = t-i + -ts 2 y 

4 2k 

Mo = Z -ts 2 - X -J z . (26) 
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FIG. 3: Critical hopping strength t c (<5,A) for quantum phase transition between Kondo singlet (t<t c ) and antiferromagnctic 
(t>t c ) phases at zero field. For (<5,A) = (0,A) or (5,1) t c is obtained from a closed analytical expression (Eq.(25) in Ref. Q) 
(full lines), for the other cases it is calculated from the zero-field version of Eq. (|23(l (dashed lines). Top: t c -dependence on 
anisotropy A of local spins . In the Ising case (A' — > 1) t c for 5=0 does not vanish because the AF order is due to the mixing 
with the doublet separated by a finite excitation gap. Bottom: t c -dependence on anisotropy 8 of interacting spins. The value 
t c (0,l)/J:r=0.7 agrees with results from MC simulatons for the genuine Kondo lattice model in 2D. [2p|. l2ll] 



Wc define a total transverse triplet amplitude i by P = \{u + d) 2 which should not be confused with the hopping 
matrix element t of the Hamiltonian. The minimization with respect to u,d leads to a relation between the two 
amplitudes from which we obtain 



f2 _ .,,2 ! + K . 2 ,2_r,72/ 1 K 



(1 + k) a M + kJ 



or equivalently 



2 _ 2? 2 _ 2Pn 2 , 2 



u=- ^; cf = - - = n 2 u. (28) 

(1 + k) 2 (1 + k) 2 y ' 

The remaining set of singlet and triplet amplitudes (s, t, to = stq) is again determined by the solution of three coupled 
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FIG. 4: In these figures three possible cases for quantum phase transitons between Kondo singlet KS(E 9 = singlet gap) and 
antiferromagnetic AF(m s = staggered moment) phases as function of increasing field strength h are shown: (i) field induced 
AF with sequence KS -AF. (ii) destruction of AF and field induced gap opening corresponding to the sequence KS-AF. 
(iii) reentrant behaviour corresponding to KS-AF-KS. Top: Singlet-triplet excitation gap E 9 and AF order parameter m s as 
function of external field for the two cases (S, A) =(0,0) (circles) and (8, A) =(1,1) (diamonds) as function of the external field. 
A subcritical scaled value t/ t c (S, A) = 0.6 is used in both cases (t c (0,0)/Ja;=0.350, t c (l,l)/J a; =0.862). Here (g s ,g T ) = (2,0) is 
used. This plot corresponds to the field induced AF case (i). Bottom: Similar plot for (5, A) = (0,1) and an above critical t/ 
t c (5, A) = 1.1 (t/Ja,=0.77) and subcritical t/ t c (<5, A) = 0.94 (t/Ja;=0.658). In the former case (ii) one has the sequence AF-KS 
of phases (A) and in the latter (iii) a reentrant situation with KS-AF-KS sequence is observed (©). 



equations: 

- 2 _ W~Y. a F a )-(u 2 + d*) 



s 



'0 



f o = ~ n , i, OT 27 T v , 77=5 = -27- to, (29) 



f 2 



2/j,+ \{2J X - J z ) + zSts 2 
-[2fi - \{2J X + J z ) + 2 7 _f + lj2 G »} + 5s 2 fl 



zt 

These amplitudes are related to the spin expectation values by 

(S 2 ) = st + ^{u 2 ~d 2 ) 

(r z ) = ~si + ^(u 2 -d 2 ) (30) 

For zero field (fo=0) the system of equations in Eq. (f!?9"|) may be shown to reduce to the one already discussed in 
Ref. In this case k—1 and u=d, i.e., circular polarised triplet modes with equal amplitudes. 
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FIG. 5: t-h Phase diagram (critical t c (h) curves) for various cases of anisotropies (8, A) and g-factors (g s ,g T ). Here circles (o) 
correspond to (<S,A) = (0,0) stars (*) to (<5,A)=(1,1) and triangles (A) to (<5,A)=(0,1). Full symbols are obtained by m 3 =0 from 
the AF side, open symbols by E 9 =0 from the Kondo singlet side. For unequal g-factors the Heisenberg case shows saturation 
at tc/Jx — 0.4 at larger fields, in the Ising case t c continues decreasing. For the mixed case (<5,A) = (0,1) (the genuine KN 
model in 2D) and unequal g-factors the t c (h) is nonmonotonic. In each case the KS phase is below and the AF phase above the 
corresponding t c (h) curve. A decreasing t c (h) leads to the phase sequence KS-AF (field- induced AF) whereas an increasing 
t c (h) entails the opposite sequence AF-KS. In the nonmonotonic region the reentrance behaviour KS-AF-KS is observed. For 
comparison see Fig. f?] For equal g-factors in all cases t c decreases monotonically and vanishes at rather small fields. 



Finally we discuss the uniform magnetisation Mo and the staggered magnetisation M s =Mb=-M-A associated with 
AF order. By definition 

M = g s (S z ) + g T (r z ) 

M s = g s {S*)+g T {T*). (31) 

Here the staggered moment is assumed to be polarised along x. Using the expressions in Eq. (|3T)|) for (S z ), (t z ) and 
similar ones for (S x ), (t x ) we obtain with the help of Eqs. (|25I27[) 

P 

M = st [2g- + 2<?+A(^ r )] ~ 2 ff _st = 2 ff _m , 
sto 

M a = st[2g- - 2g + \( t 4)} ^ 2 5 _sf = 2g^m s . (32) 
s 

Here we used g± — \{g s ± <?t) and defined 

A = i^. (33) 

1 + K 

When both staggered and uniform component are nonzero the total magnetic moment is canted. The canting angle 
a, counted from the xy plane is then given by 

,'i \9- 

tan a 



m y -^ y+A W . 04) 



The terms ~ A are higher order corrections, e.g. ~ (h/J x ) 2 for M s . For moderate fields they may be neglected. In 
Sec. IVIII (Fig. H|) we only show the main contribution m s = si in accordance with Ref. 1. In Eq. a similar small 
field approximation may be employed for the canting angle leading to a ~ tan -1 (to /i). For equal g-factors g s =g T 
(g_=0 and ?o=0) one has a special case: M s = due to moment compensation of local and itinerant spins , although 
this is still a phase with an AF order parameter m s —si. The uniform moment Mo is now given by Mg = 2g + \P, 
oriented parallel to the external field (a = ir/2). 
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The ground state energy: condensation vs. fluctuation 

It is instructive to rewrite the ground state energy in a different way which renders a clearer interpretation of its 
individual contributions. The total mean field Hamiltonian in Eq. (|18[) may be rewritten as 

n mf = E + J2 H^( fc )K + ^]> ( 35 ) 

a=z,A,B k 

where n Q is the occupation number of bosons for each mode. The total ground state energy is then given by 

Eo = £ + ^J> Q (fc), (36) 

ka 

where E$(h) consists of three parts according to 

Eo/N = -[I(2J :c + J,)( s - 2 + i) + M (^- s - 2 )] 

+|t(/t - 1) W + {\{2J X - J z ) + fi+ Z -5ts 2 )tl (37) 

Here the first part (Eq) in Eq. (|36p is the condensate contribution. The second always positive term is due to triplet 
quantum fluctuations. According to Eq. (|37[) E is composed of three parts: The first (negative) contribution is due 
to the singlet formation. The second (positive) part is due to the field induced triplet polarisation. The last term is 
the Zeeman energy contribution. For zero field only the first term in Eq. (|3"T|) is present [![. The ground state with 
Eq < will be determined by the competition of these terms. For example when h = the quantum critical point 
where AF order appears will be determined by the balance of two terms: The negative singlet formation energy and 
the positive energy of triplet quantum fluctuations. 



VI. INFLUENCE OF HIGHER ORDER TERMS 



The transformation of the Hamiltonian in Eq. (TTJ) to bond operator variables also creates third and fourth order 
terms in the triplet operators. They have been neglected in the previous analysis based on Eqs. (|6l7p . In fact it 
was argued in |1[ that third order terms do not contribute to the ground state energy and fourth order terms are 
quantitatively negligible. The latter was also found for the related bilayer Heisenberg model [23]. Hence higher order 
terms have no influence on the zero-field quantum critical properties. In this section we investigate to what extent 
this is still justified in the presence of an external field. Indeed in mean field approximation used here the higher 
order terms lead to field induced effective bilinear triplet terms which have to be added to the genuine bilinear term 
in Eq. (fTTj) which is present already at zero field. The quantitative influence of higher order terms on t c (h) or h c (t) 
is controlled by the ratio (io/s) at the critical field h c where to(ft, c ) is the field induced triplet amplitude. If it is 
still moderate compared to the singlet amplitude s(h c ) then higher order terms have negligible influence. The field 
dependence of these amplitudes as obtained from Eqs. (|23l29p is shown in Fig. [2] for a typical case. 

The contribution of third order in t na (n= site, a=x,y,z) to the interaction term is given by [l|: 




(n,m) 

+ 5[s n t niZ + t n , z Sn){t\ri,x^"rn,y ~ $ m ,y^m,x)]- (38) 

which has to be added to Eq. ([8]) . We are only interested in the influence of this term on the quantum critical lines 
h c (t) or t c (h) to be discussed in Sect. IVlil As explained there this may be achieved both from the paramagnetic (KS) 
and AF side of the critical line. Here we consider only the former for reasons explained at the end of the section. 
In this case the circular triplet amplitudes u = d = 0. The mean field approximation to H\ ' then contains only 
terms proportional to to and no constant contribution to Eo appears. Transforming these terms to circular triplet 
coordinates one finally obtains another bilinear contribution with 

n ?lf = -^ zS si J2{uiu n - did n ) + ^si a ( u l u ™ ~ d l d m)- (39) 

n ~ (rim) 
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This contribution is field induced since for small fields io(h) ~ h. The first single site term in Eq. l|39|) has the same 
structure as the 7+ part of the Zeeman term and may be accommodated by a simple (nonlinear) rescaling of the 
applied field such that 7+ is replaced by 7+ according to 

%(h) = - 9+ hUh); /.(/o = i + f*t^j^ (40) 

Here f s (h) is the rescaling function for the applied field. Note that for 5 = one has f s (h) = 1, i.e. no rescaling will 
occur in this case. The second contribution in Eq. (|39p is an interaction to which the same (transverse) Bogoliubov 
transformation as before may be applied. These terms then have a simple effect: In Eq. (1171) one has to replace 
d u ,d(k) — > d u ^d{k) ± (t/2)si j(k). However only the average d±(k) = (d u (k) + dd(k))/2 of the auxiliary functions, 
which is unchanged, enters the expression for u>A,B(k) in Eq. (|19p . Therefore u>A,B(k) is also unchanged by the 
third order contribution, As a result the total energy Eo in Eq. (|2"Tj) will be exactly the same as before. Also the 
selfconsistent equations Eq. (|2"3")l will be unchanged. Note that in these equations 7_ is not rescaled. The only effect 
of the third order terms is therefore the above rescaling of the external field in the 7+ Zeeman term leading to the 
modified transverse mode frequencies 

Q A>B (k) =o; A ,fl(A)±7 + . (41) 

Therefore the critical field h c of the quantum phase transition which is defined as the field where one of the above 
modes (f2+(fc) for h > 0) vanishes at the AF wave vector Q=(tt, tt) will be changed by the scaling factor f s (h c ). If h" 
is the critical field without third order terms then approximately h c = h° c / f s (hP c ) is the critical field with the effect of 
third order terms included. This is only a quantitative modification which we will discuss in Sec. IVIll in connection 
with Fig.O The qualitative topology of the phase diagram will not be changed by the third order term. We note again 
that in the special case (A, 6)=(A, 0) the third order term does not have any effect at all because f s = 1. Furthermore 
for all cases g s =g T there is no field induced triplet amplitude, i.e. to = according to Eq. (|23[) leading again to f s = 1. 
We conclude that the influence of third order terms is not important for the field induced quantum critical behaviour . 

Now we discuss the effect of terms which are of fourth order in the triplet operators. They are given by [l[ 

%i ' = — {(tfiytnz ~ ^nz^ny){tl ny t mz — t\ nz t m y) + {t\ x t nz — tl lz t nx )(tl nx t mz — i^fmx) 

(n,m) 

+ ^{Aix^ny ~~ ttiytnx)(tmxt'my ~ ^tny^rnx)j • (42) 



The mean field approximation to for the nonmagnetic case (u = d = ) does not produce a constant term but 
leads to an effective bilinear contribution in circular triplet coordinates 

Wfij = l u i Um + d i dm + \{ u nd n + d n u m + u) n d) m + 4«L)]- (43) 

(nm) 

Again this is a field- induced bilinear contribution with to(h) ~ h 2 for small fields. Due to u, d — it has no longitudinal 
part. It may be diagonalised by the same transverse Bogoliubov transformation in Eq. (|14|) as before. The new mode 
frequencies including fourth order contributions are simply obtained by the replacement d± = d± + (t/2)t^j(k) and 
f± = f± — (t/A)t^j(k) in Eq. p?|) . This will lead to the modified transverse frequency 

^W 2 =^(fc) 2 [l+?oT%L (44) 

and new mode frequencies flA,B(k) = u>A,B{k) ± 7+ At the quantum critical point where £l A (Q) = this may be 
approximately written as 

Q A ,B(k) - "A,B(k)[l + \(j)h(k)], (45) 

where the index c denotes the values of singlet and triplet amplitudes at the QCP. From Fig. [2] we can estimate that 
the prefactor in Eq. (|45]) is of the order j(y) — 10~ 2 since at the QCP the triplet amplitude is still quite small 
and the singlet amplitude basically unchanged from the zero-field value. Therefore the fourth order terms lead to 
corrections of the order one per cent in the frequencies tOA,B(k) which determine the last term in the ground state 
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energy of Eq. (|2Tj) . Indeed the correction to Eq is even smaller since the momentum summation over wa bQs) in the 
last term of Eq. pTj) leads to a large amount of cancellation because j(k) is positive in one half of the Brillouin zone 
and negative in the other half. Therefore the selfconsistent equations for 5 and to are nearly unchanged and we can 
conclude that the field induced quantum critical behaviour is not influenced by the inclusion of fourth order terms. 
Our analysis shows that the quantum critical lines h c (t) are only marginally influenced by the presence of higher order 
terms in one example and have no effect in other cases. However for fields with h 3> h c one would expect that the 
higher order terms should lead to non-negligible corrections, at least in the case of third order contributions. From 
Fig. [5] one can see that for h/h c > 2 the ratio (io/s) approaches 0.5 and the bilinear approximation might become 
inadequate. In fact, by its very construction the bond-operator method is a strong coupling theory which assumes 
the dominance of the singlet state. It cannot be expected to be quantitatively correct for high fields when saturation 
of moments is approached, i.e. when singlet and triplet amplitudes become equal. In this limit it is more appropri- 
ate to start from the polarised canted state and perform a conventional spin wave expansion for the two types of spins. 

As a consequence of the above analysis we will only consider the genuine bilinear Hamiltonian Eq. (jlip as used in 
the previous sections for the following numerical calculations. This restriction has an additional reason: Although the 
critical h c (t) or t c (h) may be calculated by searching for the vanishing of the excitation gap E g from the nonmagnetic 
singlet side, it is more satisfactory to search also from the AF magnetic side for the vanishing of the staggered 
magnetisation m s and confirm the agreement. This cannot be done easily when third and fourth order terms are 
included. On the AF side the latter lead to a mixing of the longitudinal and transverse modes and a closed analytic 
form of the Bogoliubov transformation cannot be found. Therefore it is better to include only the bilinear terms and 
be aware of the trivial effect of field rescaling by the third order terms and the tiny effect of fourth order terms. This 
standpoint will be adopted for the following numerical analysis. 

VII. NUMERICAL SOLUTION FOR EXCITATION GAP AND AF ORDER PARAMETER, CRITICAL 

FIELD CURVES AND THE H-T PHASE DIAGRAM 

We shall now discuss the numerical solutions of the selfconsistent Eqs. (|23I29[) which describes the AF-KS quantum 
phase transitions, i.e., the field dependence of the Kondo singlet gap E g and the field dependence of the staggered 
magnetisation m s . The former was defined below Eq. (fTi?)) and the latter in Eq. (f3"2"|) . We shall only discuss the results 
for D=2 and 5 < 1, i.e., the easy xy-plane situation with m s _L c. In this case £lA,B(k) < Vl z (k). As shown before the 
qualitative behaviour of D=3 is similar to D=2 [l| . The main purpose is to study the dependence of the quantum 
critical point t c (h) on the external field, or cquivalcntly the dependence of the quantum critical field h c (t) on the 
hopping strength t. To check consistency the QCP has been obtained both by variation of t and of h and both from 
the paramagnetic (E g =0) as well as the antiferromagnetic (m s =0) side of the QCP. As mentioned in the previous 
section this requires restriction to the genuine bilinear Hamiltonian terms. 

Before discussing the field dependent results we show the behaviour of the quantum critical t c (A, 6) as function 
of the anisotropics in the zero field case as a starting point of our analysis. The calculation of i c (A, S) has been 
described in Ref. [l|. Here we present the results for the full range of anisotropies (Fig. [3]) from the U(l) -xy case via 
Heisenberg- SU(2) to Ising-Z2 symmetry. Note the two different scales on the left (J x ) and right (J z ) half of each 
part of Fig. [3] The quantum critical t c for the KS-AF transition reaches a singular maximum for the Heisenberg 
case and generally vanishes for the Ising case, except for 6=0. This may be understood by comparing with Fig. [T] 
For A — + oo (i.e. A' — > 1) the two singlets with S* =0 become degenerate and an arbitrary small interaction t leads 

to the AF state. However for <5=0 the S* =0 states are not connected by a matrix element of Ti.^ and hence AF can 
occur only via matrix elements between the singlet ground state and the excited doublet which requires a finite t c /3 x 
= 0.35 for AF order to occur. 

The effect of the external field for (g s ,gv) = (2,0) is shown in Fig. 0] (top) for the two cases (A, 5) = (0,0) (o) and 
(1,1) ({>)■ Their corresponding t c values may be taken from Fig. [3J In Fig. 0] (top) we choose a subcritical value 
t/t c (A, S) = 0.6 meaning a nonmagnetic singlet phase with finite E g exists for zero field. When the field is increased 
the gap is gradually reduced until it is closed at t c (h c ). For the (normalised) t value of 0.6 obviously the difference in 
critical fields h c is not significant. However, the qualitative behaviour of the spin gap E g in the two cases for h<h c 
is different with a much smaller slope of E g in the xy-case as compared to the Heisenberg case. In both cases for 
h>h c a field induced easy plane AF order parameter appears, i.e., a phase sequence KS-AF for increasing field. The 
opposite behaviour is observed in the genuine KN case (A,<5) = (1,0) with (g s ,gi-) = (2,0) in Fig. Q] (bottom). For an 
above- critical value t/t c (0, 1) = 1.1 (t/J K =0.77) the AF order parameter is gradually suppressed until it vanishes at 
hc/Jz =1.23 and the Kondo singlet phase appears, i.e., the opposite phase sequence AF-KS is realised. Finally, for 
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FIG. 6: Canting angle a of total moment moment (counted from the xy-plane) in the AF phase of two different cases (cf. 
Fig. [4| with (g a ,g T ) = (2,0). Full curve is obtained from the exact expression Eq. (|34[) and dashed curve from the approximate 
expression a ~ tan - (to/i). At the quantum critical field h c the staggered moment vanishes and the total moment is aligned 
with the external field (a = tt/2). Here h c (0,l)/J a: = 1.23 and h c (l,l)/J a; =0.55. 



slightly subcritical t/t c (0, 1) = 0.94 (t/J^O.658) an interesting reentrance sequence KS-AF-KS of phases as function 
of field is observed. It exists only in a narrow range of subcritical values 0.91 < t/t c < 1.0. As discussed in Sect. IVII 
this behaviour is robust because third order triplet contributions are exactly zero for 5=0. 

We may collect the data of quantum critical fields h c from numerous calculations such as presented in Fig. [4] and 
similar ones for quantum critical t c at fixed h in a h-t phase diagram. It is shown in the form of t c (h) curves in Fig. [5] 
for two choices of the g-factors. For (g s ,g T ) = (2,0) and the cases (A, S)= (0,0), (1,1) the scaled t c (h) is identical 
close to the zero-field QCP, i.e., the slope does not depend on the anisotropy (A,<5). For larger fields one observes 
the plateau formation in the case (A,<5)=(1,1) and the further monotonic decrease for (A, 6)= (0,0). This behaviour 
may qualitatively be understood from the field dependence of the local energy levels of 7ij shown in Fig. [TJ For 
increasing h the gap between the ground state and first excited state decreases and hence a smaller t c is necessary 
to achieve the softening of the triplet excitation at Q. For A=l this effect eventually levels off because the splitting 
of the two lowest local levels becomes constant at large field and approaches 2A3 X so that t c (h) reaches a plateau at 
large fields for nonzero A. On the other hand for A=0 the two lowest levels become asymptotically degenerate and 
hence t c (h) should approach zero for large fields. As already noted in the discussion of Fig. Q] the intermediate case 
(A, S ) = (1,0) and (g s ,g T ) = (2,0) behaves quite differently. After an initial but much weaker decrease of t c (h) it 
reaches a minimum at (h/J x ,t/3 x )=(0. 275, 0.638) and then starts to increase again. The region of the nonmonotonic 
behaviour with slightly subcritical t/t c (0,l) corresponds to the region where reentrance behaviour KS-AF-KS is 
observed in Fig. 2] For above critical t/t c (0,l), once t c (h)>t the initial AF phase is suppressed and we obtain the 
sequence AF-KS, i.e. a field induced KS phase corresponding to Fig.. [4] (bottom). In all three cases with (g s ,g r ) = 
(2,0) one has a finite t c (h) for finite fields. 

This is different if we chose equal g-factors (g s ,g T ) = (2,2). For this choice the lowest two local levels cross at 
h=0.25 (A=0) and at h=0.5 (A=l) (not shown in Fig. [I]). This is due to the fact that S' commutes with the 
Hamiltonian and hence there is no level repulsion of the two S*=0 states. Consequently at this point the quantum 
critical t c (h) vanishes in all three cases of (A, (5). It decreases essentially linearly from its maximum value at h=0. The 
phase diagram in Fig. O shows a collection of the critical t c (h) curves for the six cases considered. The corresponding 
gapped KS and AF ordered phases are below and above each curve respectively. We comment again on the influence 
of third order terms discussed in Sect. IVII on the phase diagram of Fig. [5] Remembering that all cases with 8 = or 
g s = g T are unaffected by the third order contributions, only one of the critical field curves in Fig.[5](* with (g s ,g r ) = 
(2,0)) may be influenced. In this case third order terms lead to an almost field independent and moderate rescaling 
of h c by a factor ~1.2 according to Eq. (|40|) and Fig. [2] 

Finally in Fig. [S] we show the dependence of the tilting angle of the total moment out of the xy-plane for two 
possible cases with AF order. It is calculated from Eq. (fM)) (full line) and a simplified (low field) expression given in 
the caption. Close to the quantum critical field the staggered moment vanishes and the total moment is then alligned 
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with the external field (|| z), meaning a — k/2. 

VIII. DISCUSSION 

As we have mentioned in the introduction the interest for the Kondo necklace (KN) model stems from two 
facts: (i) it is thought to be a simplified version of the original Kondo lattice (KL) model with itinerant electrons 
and (ii) it is expected to describe some basic features of the observed quantum phase transitions in real Kondo 
materials. Although there is no mapping or exact equivalence of the two models the agreement of their qualitative 
features by simply identifying the intersite interaction (KN) and hopping (KL) energies (both given by t) is generally 
quite good. As discussed in [3, [TtJ the ID Kondo necklace model does not exhibit quantum critical behaviour, 
this agrees with earlier numerical work on the KL model [23[ where the disordered spin gapped state always prevails [1[ • 

We discuss the more interesting 2D case where quantum phase transitions in the KN model are possible [l|, [ijj • 
One can compare the values of the quantum critical t c /J (we set J= in this section to comply with literature 
conventions) obtained for the KN model from Fig. [3] with those for the KL model obtained in MC simulations 
obtained previously 0, [U H3] . In the spirit of Doniach's replacement in ID we should compare with the KN model 
for xy-type interactions, i.e. with the (6, A) = (0, 1) case in Fig. [3] which has a t c /J=0.7. The MC simulations 
for the 2D SU(2) fermionic KL model lead to quantum critical parameters (average values) t c /J= 0.71 [2fJ , 0.69 
[2ll l2~ij and 0.68 [lfij], very close to the appropriate value in Fig. [3J (upper panel: center of full line; lower panel: 
left corner of full line) . Note that the critical t c 's corresponding to full lines in Fig. [3J are obtained from a simple 
analytical formula within the bond-operator mean field approach [l|. Further support for this scenario comes from 
exact diagonalisation (ED) results for the KL model on small clusters [25[ using open boundary conditions. The total 
average local moment (fJ.f oc ) = {{ti + Si) 2 ) calculated as function of J/t shows a reduction to half of the free moment 
size due to Kondo singlet formation in the lattice for the value t c /J = 0.67. If it is interpreted as the critical value in 
the thermodynamic limit this is again close to the above results. This near equality of analytical KN and numerical 
KL results seems to suggest that even in 2D the SU(2) KL model is better described by the U(l) xy-type KN model 
with (6, A) = (0,1) since the SU(2) KN model with (5, A)=(l,l) has a considerably larger critical value t c /J = 0.862. 

The fermionic KL model is obtained from the periodic Anderson model (PAM) by a Schrieffer- Wolff trans- 
formation which eliminates charge fluctuations of f-electrons. The resulting Kondo coupling J is then given in 
terms of the Anderson model parameters by J = —UV/ef(et + U) where e/ < is the f-level position with 
respect to the Fermi level at zero energy, U is the on-site f-electron repulsion and V their hybridisation strength 
with conduction electrons. The periodic Anderson model may be studied with the dynamical mean field theory 
(DMFT) 26]. This method has been used more recently to investigate its quantum critical properties [27j which 
should be related to those of the Kondo lattice and hence also KN models. The method is formally a D=oo 
approximation but has nevertheless been used for studying electron correlations in finite dimensional lattices 
such as the 3D PAM in [27|. Using the value of J from the Schrieffer- Wolff transformation the extrapolated T= 
phase diagram suggests that the critical hopping strength at the QCP is given by t c /J ~ 0.43. This is again 
in reasonable agreement with the mean-field bond operator result for the 3D KN model which predicts t c / J — 0.38 [l[. 

While the quantum critical behaviour as function of control parameter t/J is well presented, there are few 
microscopic investigations concerning the field induced quantum critical behaviour of KL/KN type models such as 
provided here. The fermionic 2D KL model in an external field with equal g-factors corresponding to (g s ,g T ) — (2,2) 
was studied in 

RefQJ by using variational and MC methods. The numerical results for J/t =3 suggest that the 
transition between Kondo and AF ordered regime takes place in a field region (B= 2h in notation of Ref. ITil ) between 
B~/t = 1.0 and n B+/t =2.25. In this regime both Kondo-like features and transverse magnetic order coexist. From 
the present calculations of the 2D KN model again with (S, A)=(0,1), g-factors as above and using the data shown in 
Figs. l3l5l one obtains B c /t = 1.86 (B c /J=0.62) which is close to the average 1.63 of the above B^ jt values. This value 
for the critical field was also confirmed by exact diagonalisation (ED) results for the KL model on small clusters [25[ 
where the breaking up of on-site singlet formation was observed above the value B c /J ~ 0.5. We conclude that KL 
and KN models also seem to exhibit similar field induced quantum critical behaviour. Another matter is the scaling 
behaviour of singlet gap and staggered magnetisation around the critical field. The present mean field type theory 
does not give sufficient insight into that issue. It has sofar been treated within continuum field theories involving 
order parameter fluctuations around the critical field [28j |. 

The present bond operator mean field results for the KN model are apparently consistent with known numerical 
results for the fermionic Kondo lattice models. It is less clear whether they are useful for the interpretation of 
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experimental results, especially for the field induced quantum phase transitions. Since charge degrees of freedoms 
are eliminated in the KN model it is strictly speaking more appropriate for the Kondo-insulator compounds which 
have a hybridisation gap due to half filled conduction bands. However we ignore this subtlety in the following and 
also apply it to the magnetism of metallic Kondo compounds. The zero field quantum critical behaviour as function 
of the control parameter t/J can experimentally be mimicked by applying hydrostatic or chemical pressure (by 
substitution of elements). This changes mostly the hybridisation and hence the Kondo coupling J. In the accessible 
pressure regime one may assume that t/J varies linearly with pressure. In this manner it is feasible to drive AF 
ordered heavy fermion systems to the quantum critical point where they become nonmagnetic heavy Fermi liquids. 
Since hydrostatic (positive) pressure generally increases J it tends to suppress the AF phase while with suitable 
substitution of elements (negative) chemical pressure may decrease J which favors AF order. Thus the AF QCP may 
be approached from both sides. There are many examples to be found especially among Ce- compounds, for a review 
see In most cases of Ce-compounds, however, the AF QCP may not be reached directly because it is enveloped 
by a 'dome' of the superconducting phase and the critical pressure has to be obtained from extrapolation to T=0. 
While the KL/KN type models discussed above naturally suggest the AF QCP there is no direct way of experimental 
determination of the critical (t/J) c since its relation to the experimentally accessible critical pressure p c (or critical 
concentration of substituent) is unknown. In fact most of the interest on pressure induced QCP's focuses more on the 
scaling exponents of various quantities with respect to the distance (p-p c ) to the QCP. This is extensively discussed 
in the reviews cited in the introduction. 

For the field induced quantum phase transition much less experimental results are available. A review of 
materials investigated is given in Ref. [3(| A recent example of a field induced destruction of AF order is the 
tetragonal YbRli2Si2 compound [3l[. At ambient pressure and zero field it has an AF order which is indeed of 
the easy plane (xy)-interaction type (5=0). The AF order is destroyed at 0.66T for field along the hard c-axis 
and a nonmagnetic Fermi liquid state appears above the critical field . If one assumes that the local Kondo 
interaction which results from the Schrieffer- Wolff mechanism is more of the isotropic nature (A=l) it qualitatively 
corresponds to the AF-KS scenario shown in the lower panel of Fig. 2] As mentioned before the reentrant scenario 
shown also in this figure exists only in a narrow parameter range and it is perhaps understandable that no such 
example is known. What is really surprising is the following observation: Most of the parameter cases in Fig. [5] 
would predict the KS-AF sequence, i.e. the field-induced AF order out of the Kondo phase. This is a natural 
consequence that in most cases a field reduces the singlet-triplet gap and supports the onset of magnetic order 
as shown in the upper panel of Fig. We would like to stress that this is also the only phase sequence found 
in the fermionic KL model [lij]. However field induced AF magnetic order is not easily found in heay fermion 
compounds. For example CeNi2Ge2 at ambient pressure and zero field is in the Kondo singlet phase and rather 
close to the AF quantum phase transition which may be achieved by appropriate substitution of Ni. Thus it is a 
complementary case to YbRli2Si2. However in an external field there is no field induced transition to transverse 
AF order as one might naively expect from the model discussion, rather it is driven further away from the AF 
QCP as specific heat and resistivity measurements suggest |3lj . There is however one known example of a heavy 
fermion system which exhibits field induced AF order, though complicated by the appearance of superconductivity. 
The compound CeRhlns is already an AF at ambient pressure [32[. Application of hydrostatic pressure moves it 
to the nonmagnetic (KS) side where it also becomes superconducting. Additional application of a magnetic field 
destroys the superconductivity and leads to the reestablishment of a field induced AF order in a wide range of 
pressure form 1.4 to 2.4 GPa [32|]. This is indeed the KS-AF sequence of phases corresponding to Fig. U (top) which 
appears in the majority of cases studied here. Experimentally, however, it seems to be the one most rarely encountered. 

These examples suggest that while the Kondo lattice or necklace type models have some qualitative and instructive 
properties which are relevant for magnetic quantum phase transitions they may not be able to explain some experi- 
mental observations in heavy fermion compounds. It is well possible that the complete elimination of charge degrees 
of freedom is too radical by suppressing hybridisation effects of moments with conduction electrons and thus favoring 
the magnetic phases. Therefore the more fundamental PAM model, including appropriate local and conduction elec- 
tron Zeeman terms may be a more promising starting point for future investigations of field induced quantum phase 
transitions in heavy fermion compounds. 

IX. CONCLUSION AND OUTLOOK 

In this work we have investigated the general anisotropic Kondo necklace model and its quantum phase transi- 
tion from paramagnetic singlet phase to AF phase and vice versa in an external field. We have derived a bosonic 
Hamiltonian using the bond operator representation of local and interacting spins and diagonalised the bilinear part 
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by a generalised Bogoliubov transformation to circular polarised modes to include the effect of the external field. 
Higher order parts of the Hamiltonian were found to be insignificant for the quantum critical behaviour. In mean 
field approximation one obtains a total ground state energy that depends on the singlet amplitude s, staggered and 
uniform triplet amplitudes t, to and the chemical potential. Minimization leads to two sets of selfconsistent equations 
for Kondo singlet and AF phase, respectively. 

From their numerical solution we find that in most investigated cases the t c value of the quantum critical point 
decreases monotonously with field strength leading to the KN-AF phase sequence for increasing field. In the case 
of unequal g-factors and (A, S) = (1,1) it levels off at a plateau value t c (h)/t c (A, S) ~ 0.4 while for (A, S) = (0,0) 
it continues to decrease at larger field. However for the genuine KN case with (A, 6) = (1,0) the t c (h) curve is 
nonmonotonic with a minimum at h/ J x ~ 0.275 at small fields and continues to increase for larger fields. Depending 
on the size of t this implies two possibilities: For above critical t a suppression of AF and the opening of a spin 
excitation gap for h > h c (AF-KN sequence) or the reentrance behaviour (KN-AF-KN sequence) of the gapped singlet 
phase for slightly subcritical t (as compared to t c (0)). For equal g-factors a crossing of local levels occurs due to a 
conserved S' and the critical t c reaches zero already at relatively small fields for all investigated cases of (S, A). Since 
our treatment is of mean-field type we will get mean-field exponents for the spin gap E g and magnetisation m s at 
the field-induced QCP. This is apparent from FigUand has not been discussed further. Improvement for the critical 
exponents will require a selfconsistent renormalisation theory for the triplet excitations close to the quantum critical 
point. 
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